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A  3D  CFD  model  for  two-stage  entrained  flow  dry  feed  coal  gasifier  with  multi  opposite  burners  (MOB) 
has  been  developed  in  this  paper.  At  each  stage  two  opposite  nozzles  are  impinging  whereas  the  two 
other  opposite  nozzles  are  slightly  tangential.  Various  numerical  simulations  were  carried  out  in  standard 
CFD  software  to  investigate  the  impacts  of  coal  and  oxidant  distributions  between  the  two  stages  of  the 
gasifier.  Chemical  process  was  described  by  Finite  Rate/Eddy  Dissipation  model.  Heterogeneous  and 
homogeneous  reactions  were  defined  using  the  published  kinetic  data  and  realizable  k-s  turbulent  model 
was  used  to  solve  the  turbulence  equations.  Gas-solid  interaction  was  defined  by  Euler-Lagrangian  frame 
work.  Different  reaction  mechanism  were  investigated  first  for  the  validation  of  the  model  from  pub¬ 
lished  experimental  results.  Then  further  investigations  were  made  through  the  validated  model  for 
important  parameters  like  species  concentrations  in  syngas,  char  conversion,  maximum  inside  tempera¬ 
ture  and  syngas  exit  temperature.  The  analysis  of  the  results  from  various  simulated  cases  shows  that 
coal/oxidant  distribution  between  the  stages  has  great  influence  on  the  overall  performance  of  gasifier. 
The  maximum  char  conversion  was  found  99.79%  with  coal  60%  and  oxygen  50%  of  upper  level  of  injec¬ 
tion.  The  minimum  char  conversion  was  observed  95.45%  at  30%  coal  with  40%  oxygen  at  same  level.  In 
general  with  oxygen  and  coal  above  or  equal  to  50%  of  total  at  upper  injection  level  has  shown  an  opti¬ 
mized  performance. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Depletion  of  conventional  energy  resources  like  fossil  fuels  have 
diverted  the  focus  of  engineers  and  researchers,  around  the  globe, 
to  develop  the  energy  efficient  systems  particularly  for  industries. 
Advancements  of  new  economic  and  efficient  technologies  like  co¬ 
generation  [1  ],  combined  cycles  [2,3]  and  gasification  [4]  are  one  of 
the  prime  concern,  nowadays  for  developing  countries  to  meet  the 
socio-economic  sustainability.  Coal  gasification  is  a  complex  pro¬ 
cess  in  which  fuel  is  converted  to  combustible  gases  (CO  and  H2) 
in  the  presence  of  water  and  oxygen,  at  elevated  pressures  and 
temperatures.  Product  gas  is  known  as  “syngas”  which  can  be  used 
as  town  gas,  to  make  electric  power,  methanol,  liquid  fuels,  urea, 
ammonia,  etc.  [5  .  Coal  gasification  is  highly  environment  friendly 
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as  it  has  the  capability  to  achieve  extremely  low  SOx,  NOx  and  par¬ 
ticulate  emissions  [6  . 

For  a  successful  and  efficient  coal  gasification,  a  deep  under¬ 
standing  of  all  physical  and  chemical  changes  involved  in  the  pro¬ 
cess,  is  required.  Numerical  simulation  provides  an  easy,  cost 
effective  and  reliable  way,  to  study  the  effect  of  various  controlling 
parameters,  like  coal  composition,  oxidant  to  fuel  ratio,  residence 
time  of  fuel  particles  in  the  system,  temperature  and  pressure  of 
the  gasifier,  rate  of  chemical  reactions,  etc.  on  the  gasification  pro¬ 
cess.  For  last  20  years,  the  CFD  modeling  remained  an  effective  tool 
to  simulate  and  visualize  the  gasification  process  [7-12  .  Different 
mathematical  models  have  been  developed  for  various  sub¬ 
processes  involved  in  the  coal  gasification.  For  instance  some 
researchers  have  developed  Euler-Euler  model  to  define  the  gas 
and  solid  flows  [6,13-16]  whereas,  the  others  have  used  Euler- 
Lagrangian  frame  work  to  setup  the  gas  and  solid  interactions 
[17-20].  Probability  Density  Functions  (PDF)  [10,21,22]  model 
and  Finite  Rate  Chemistry/Eddy  Dissipation  models  [14,23-25] 
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Nomenclature 

P*  Pp 

density,  density  of  particles  (kg/m3) 

Qr 

heat  flux  for  radiation  heat  (J/m2  s) 

u,  up 

velocity,  velocity  of  particles  (m/s) 

G 

incident  radiation 

Sm,  Sj , 

Sh,  Sr  source  terms  for  mass,  momentum,  energy  and 

C 

coefficient  of  function  for  linear-anisotropic  phase 

species 

a 

absorption  coefficient 

cp 

specific  heat  at  constant  pressure  (J/kg  K) 

(7s 

scattering  coefficient  (m_1) 

Tij 

symmetric  stress  tensor 

G 

Stefan-Boltzmann  constant 

T 

temperature  (I<) 

emissivity 

A 

turbulent  thermal  conductivity  (W/m  K) 

Wj,r 

net  production  rate  of  species  i  through  chemical 

Cj 

mole  fraction  of  specie  j 

reaction  (K  mol/m3  s) 

Di 

diffusivity  (m2/s) 

Mj 

molecular  weight  of  specie  j 

Fd 

drag  force  (kg  m/s) 

Vlr 

Stoichiometric  coefficient  for  product  j  in  reaction  r 

Fx 

force  vector  component  along  x-axis  (kg  m/s2) 

if. 
i  r 

Stoichiometric  coefficient  for  reactant  i  in  reaction  r 

P 

dynamic  viscosity  (N  s/m2) 

[C] 

molar  concentration  of  specie  (I<  mol/m3) 

k 

kinetic  energy  for  turbulence  (m2/s2) 

rj" 

rate  exponent  for  product  specie 

8 

dissipation  rate  of  turbulence  (m2/s3) 

rj' 

rate  exponent  for  product  reactant  specie 

Pt 

turbulence  viscosity 

kf 

forward  reaction  rate  constant 

c. 

viscosity  constant 

A 

pre-exponential  factor  (consistent  units) 

Gk 

mean  velocity  gradients 

B 

temperature  constant  (dimensionless) 

Dt 

diffusion  coefficient  for  turbulence  (m2/s) 

E„ 

activation  energy  for  reaction  (J/I<  mol) 

Prt 

Prandtl  number  for  turbulence 

R 

universal  gas  constant  (J/I<  mol  I<) 

Sct 

Schmidt  number  for  turbulence 

are  the  most  common  methods,  used  so  far,  to  define  the  chemistry 
of  the  gasification  reactions. 

The  effects  of  major  design  and  operating  parameters  like  particle 
density/size,  oxygen/coal  ratio,  injector  nozzle  design,  number  of 
nozzles,  gasification  reactions,  etc.,  have  been  investigated  by  vari¬ 
ous  researchers  through  CFD  simulations  [26-31  .  Entrained  flow 
gasifiers  (EFSG)  are  considered  to  be  the  best  in  advanced  integrated 
gasification  combined  cycle  (IGCC)  power  plants,  due  to  high  carbon 
conversion  efficiencies  with  short  residence  time  [32].  Opposite 
Multi  Burners  (OMB)  in  entrained  flow  gasifiers  are  grown  popular 
nowadays  for  the  better  fuel  conversion  due  to  high  collision  rates 
for  solid  particles  and  resolving  the  issues  of  deposition  of  slag/ash 
on  gasifier  walls  [33-35  .  In  the  OMB  there  are  two  sets  of  nozzles 
exactly  opposite  to  each  other  which  produce  impinging  flows.  Tang 
et  al.  [36]  has  partially  modified  the  design  of  OMB  by  changing  one 
set  of  impinging  nozzles  into  tangentially  opposite  to  each  other  and 
he  has  introduced  a  second  stage  to  develop  flameless  oxidation 
conditions.  Two-stage  entrained  flow  gasifier  with  impinging  and 
tangential  nozzles  produces  better  results  as  compared  to  OMB  gas¬ 
ifier  with  all  impinging  nozzles,  but  no  study  has  been  conducted,  so 
far,  to  investigate  the  effects  of  fuel/oxidation  amounts  at  each  level 
on  the  overall  performance  of  the  gasifier. 

The  objective  of  present  work  is  to  study  the  effects  of  coal  and 
oxidant  distribution  effects  between  two  stages  on  important 
parameters  of  gasifier  (like  temperature  and  syngas  composition, 
etc.)  through  numerical  simulations.  A  3D  CFD  model  of  two-stage 
coal  gasifier  with  impinging  and  tangential  nozzles  is  developed  in 
which  Euler-Lagrangian  frame  work  and  Finite  Rate/Eddy 
Dissipation  model  are  used.  In  the  first  part,  few  reaction  mecha¬ 
nisms  have  been  investigated  to  validate  the  model  from  experi¬ 
mental  results.  After  the  validation  of  the  model  the  effects  of 
coal  and  oxygen  distribution  have  investigated  with  fixed  oxy¬ 
gen/coal  ratio  on  the  performance  of  lignite  coal  gasifier. 


2.  Numerical  method 

2.2.  Description  of  physical  system 

An  oxygen-blown,  entrained  flow  coal  gasifier  for  Chinese  high- 
ash  lignite  coal  is  used  in  this  study  as  shown  Fig.  1.  The  proximate 


and  ultimate  analysis  of  coal  is  listed  in  Table  1.  The  particle  diam¬ 
eters  are  fitted  to  the  Rosin-Rammler  distribution  (minimum  dia  is 
4  pm,  maximum  dia  is  125  pm,  and  mean  dia  45.6  pm). 

The  inner  radius  of  experimental  gasifier  is  0.32  m  and  its 
height  is  0.55  m.  The  coal  feeding  rate  is  taken  as  7.2  kg/h  [37]. 
Four  nozzles  are  installed  at  two  levels  of  gasifier  chamber  named 
as  AA'  Level  (Top)  and  BB'  Level  (Bottom).  Two  opposite  nozzles, 
exactly  at  the  central  axis  line  for  each  levels,  are  impinging  in  nat¬ 
ure  whereas  the  other  two  nozzles  are  slightly  at  a  side  of  the  axis 
and  produce  a  swirl  in  the  flow  and  called  tangential  nozzles.  At 
top  (AA'  Level)  the  coal  is  injected  with  carrier  inert  gas  (N2)  from 
opposite  impinging  nozzles  whereas  pure  oxygen  is  injected 
through  tangential  nozzles.  However  at  bottom  (BB'  Level)  the  coal 
is  injected  with  02  from  impinging  nozzles  (Fig.  1).  The  amount  of 
oxygen  is  calculated  for  all  the  cases  on  the  basis  of  oxygen  to  coal 
(O/C)  ratio  with  fixed  coal  feeding.  The  amount  of  coal  and  oxygen 
feedings  at  individual  levels  are  calculated  on  the  basis  of  selected 
coal  or  oxygen  distribution  at  each  level  respectively. 

2.2.  Computational  models 

In  present  work,  numerical  study  is  carried  out  with  3D,  steady 
and  incompressible  turbulence  flow  with  heterogeneous  and 
homogenous  reactions.  Therefore,  time-averaged-steady-state 
Navier-Stokes,  mass  momentum  and  energy  and  species  equations 
have  solved.  The  governing  equations  for  the  system  are  as  follows 
[5]. 


( PUij )  —  $m 

(1) 

dP  d 

(pm)  -  psj  dx.  +  dx.  (t,j  pm]) + Sj 

(2) 

( PCpUtT )  -  ^  f  A  gxi  pCpU'jT')  +  p<P  +  Sh 

(3) 

(4) 

where  Ty  is  the  symmetric  stress  tensor  and  pujuj  is  the  Reynolds 
stress.  The  realizable  k-s  turbulence  model  is  also  employed  to 
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Table  1 

Properties  of  Chinese  coal. 


Proximate  analysis  (w/%) 

Ultimate  analysis  (war/%) 

Moisture 

1.16 

C 

63.51 

Ash 

23.10 

H 

4.19 

Volatile 

27.03 

O 

6.39 

Fixed  carbon 

48.21 

N 

1.02 

S 

0.63 

Heating  value 

32.91  MJ/kg 

Note:  values  are  taken  from  earlier  work  [37]. 


solve  the  turbulent  flow.  The  turbulence  kinematic  viscosity  is  cal¬ 
culated  by: 

Ht  =  Pcf  <fc2/e  (5) 


tydQ 

Sct  dX[ 


where  Prt  (=0.85)  is  the  turbulence  Prandtl  number  and  Sct  (=0.7)  is 
the  turbulence  Schmidt  number. 

Lagrangian  approach  has  been  used  to  calculate  the  motion  of 
particles  by  using  Discrete  Phase  Model  (DPM).  The  trajectories 
of  coal  particles  are  predicted  in  DPM  by  integrating  the  force  bal¬ 
ance  on  the  coal  particle,  when  they  move  through  the  continuous 
phase  of  the  fluid  [40].  This  force  balance  equates  the  coal  particle 
inertia  with  the  forces  acting  on  the  coal  particle,  and  can  be  writ¬ 
ten  (for  the  x  direction  in  Cartesian  coordinates)  as  [41  . 

=  Fd(U  ~  Up)  +gx^-—^-+Fx  (10) 

Lit  fjp 


where  CAi  is  the  viscosity  constant,  k  the  turbulence  kinetic  energy 
and  s  the  turbulence  dissipation  rate,  k  and  s  are  calculated  from 
the  following  transport  equations  [38  : 


d_ 

dXi 


d_ 

dXi 


dW 

dXi 


+  G/< 


p£ 


d_ 

dXi 


d 

dXi 


d£ 

dXi 


+  Cle  Gk 


k 


p2 

C2sGkJ 


where  Gk  is  the  generation  of  turbulence  kinetic  energy  due  to  mean 
velocity  gradients.  The  constants  in  the  transport  equations  are 
selected  from  Launder  and  Spalding  [39  and  they  are  given  by 
CAi  =  0.09,  Cig  =  1-44,  C2e  =  1.92,  ak  =1.0  and  cr£  =  1.3.  The  turbulence 
heat  conductivity  (2)  and  diffusion  coefficient  (D)  in  Eqs.  (3)  and  (4) 
are  given  by: 


pCpU'J' 


*dT  _  pt  dT 
dXi  ~  °p  Prf  dXi 


The  interaction  between  the  discrete  phase  and  the  continuous 
phase  is  also  taken  into  account  by  treating  the  heat  and  mass 
losses  of  the  particles  as  the  source  terms  in  the  governing 
equations. 

P-1  model  has  adopted  to  calculate  the  radiant  heat  in  the  gas¬ 
ifier  [9,10,13,41  .  In  PI  model,  the  local  radiation  intensity  is 
calculated  by  the  equation 

-Vqr  =  aG  -  4aGaT4  (11) 

where 


Qr 


1 

3(a  +  ers) 


where  a  is  the  absorption  coefficient,  as  is  the  scattering  coefficient, 
G  is  the  incident  radiation,  C  is  the  linear  anisotropic  phase  function 
coefficient,  and  a  is  the  Stefan-Boltzmann  constant. 
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2.3.  Combustion/gasification  model 

Species  transport  model  (Eq.  (4))  is  used  for  the  gasification 
reactions  chemistry.  This  modeling  approach  gives  an  option  to 
define  the  important  reactions  and  set  their  kinetic  parameters. 
During  the  high  temperature  environment  of  coal  gasification, 
the  coal  will  be  decomposed  into  volatiles,  char  and  ash  [42  .  The 
compositions  released  from  the  coal  can  be  expressed  by  the  fol¬ 
lowing  equilibrium  equation  43]. 

Coal  ->  ai  Volatiles  +  a2  H20  +  a3  Char  +  a4  Ash  (13) 

After  undergoing  fast  heating,  the  hot  flow  around  coal  particles 
will  trigger  a  number  of  physical  and  chemical  reactions  [29  .  The 
reactions  include  the  devolatilization  of  coal,  the  combustion  of 
volatiles  and  unburned  char  as  well  as  the  gasification  of  the  char 
[44  .  In  present  study,  the  volatiles  are  lumped  into  one  volatile  gas 
species  C1.37H4.58O0.44  and  are  calculated  from  proximate  and  ulti¬ 
mate  analysis  of  coal.  The  volatiles  release  are  described  by  a  two- 
step  devolatilization  model  [45]  and  it  is  given  as  follows: 

Coal^(  1  -Yi)  x  Char^Y  1  x  Volatile  (for  low  temperature)  (14) 

Coal^(  1  -Yh)  x  Charh  +  Yh  x  Volatile  (for  high  temperature)  (15) 

where  Y  is  stoichiometric  coefficient.  At  low  temperature  Eq.  (14)  is 
dominated  whereas  Eq.  (15)  shows  a  higher  rate  at  high  tempera¬ 
ture.  The  reaction  kinetic  equations  are  as  follows: 

^ =(k,Y,+khYh)Coal  (16) 

ki  =Aiexp(-Ei/RTp)  (17) 


kh=Ahexp(-Eh/RTp)  (18) 

where  V  denotes  the  mass  fraction  of  volatiles,  k  is  the  reaction  rate 
constant,  A  the  pre-exponential  factor,  TP  the  coal  particle  tempera¬ 
ture  and  E  the  activation  energy  of  the  reaction.  The  values  of  Y/,  /</, 
Yh,  kh,  E{  and  Eh  are  obtained  from  previous  studies  which  have  been 
done  on  similar  type  of  coal  [45,46]  and  they  are  listed  in  Table  2. 

When  char  is  produced  from  coal  devolatilization,  CO  and  H2 
can  be  generated  from  char  gasification.  There  are  various  reac¬ 
tions  selected  by  different  researchers  to  define  the  gasification 
reaction  mechanism  [5,6,13,14,17-20,23,25  .  In  this  paper,  preli¬ 
minary  simulations  were  carried  out  to  find  the  best  reaction 
mechanism  among  different  reaction-sets.  The  details  of  those 
cases  are  tabulated  in  Table  2. 

The  various  reaction  mechanisms  are  modeled  with  following 
chemical  species:  C(s),  02,  N2,  CO,  C02,  H20,  H2  and  volatiles.  Finite 
Rate/Eddy  Dissipation  rate  model  has  been  used  to  calculate  the 
rate  of  formation  for  each  specie  and  update  the  source  term  Sr 
in  Eq.  (4)  by  the  following  expression. 


Sr  =  (19) 

j=i 


kf  =ATBe(-E°/RT)  (21) 

In  the  above  equation,  the  forward  reaction  rate  constant  kf  is 
established  based  on  Arrhenius  law;  A  is  the  pre-exponential  fac¬ 
tor,  B  the  temperature  exponent  and  Ea  the  activation  energy  of 
reaction.  The  values  of  A,  B  and  Ea  for  various  reactions  are  obtained 
from  earlier  studies  and  given  in  fable  3. 


Table  2 

Various  preliminary  cases  to  optimize  the  reaction  mechanism. 


Reactions 

Simulation 

cases 

A 

B 

C 

D 

E 

F 

Vol  +  2.295  02  ->  1.37  C02  +  2.29  H20  (volatiles  complete  combustion) 

v* 

Vol  +  1.61  02  ->  1.37  CO  +  2.29  H20  (volatiles  partial  combustion) 

v* 

C(S)  +0.5  0; 

2  ->  CO  (char  partial  combustion) 

is* 

V* 

C(s)  +  02  - 

C02  (char  complete  combustion) 

is* 

C(s)  +  C02  - 

■»  2CO  (gasification,  Boudourad  reaction) 

is* 

is* 

is* 

V* 

C(s)  +  H20  - 

->  CO  +  H2  (gasification) 

is* 

is* 

V* 

CO  +  0.5  02 

->  C02  (CO  combustion) 

is* 

H2  +  0.5  02 

->  H20  (H2  combustion) 

is* 

is* 

u* 

v* 

CO  +  H20  <-> 

C02  +  H2  (water  gas  shift  reaction) 

is* 

v* 

Table  3 

Selected  kinetic  parameters  for  devolatilization  and  gasification/combustion  reactions. 

Devolatilization  [45,46] 

V/:  Yh 

0.3;  1 

kh  kh  (s”1) 

2  x  105;  1.3  x  107 

Er,  Eh  (kj  mol-1) 

104.6;  167.4 

Combustion/gasification  reactions  [31] 

A 

B 

Ea  G  K  mor1) 

Heterogeneous  reactions  ( solid-gas  phase) 

C(s)  +  02  -  C02 

0.002 

0 

7.9  X  107 

C(s)  +  0.502  ->  CO 

0.052 

0 

6.1  X  107 

C(s)  +  C02  ->  2CO 

242 

0 

2.75  x  10s 

C(s)  +  H20  ->  co  +  h2 

426 

0 

3.16  x  108 

Homogeneous  reaction  (gas  phase) 

CO  +  0.502  ->  C02 

2.239  x  1012 

0 

1.7  x  10s 

H2  +  0.5  02  ->  H20 

6.8  x  1015 

0 

1.68  x  10s 

CO  +  H20  ^  C02  +  H2  (WGS  reaction)  / 

2.75  x  1010 

0 

8.38  x  107 

b 

2.65  x  10-2 

0 

3.96  x  103 

C1.37H4.5gO0.44  (volatile)  +  1.61  02  ->  1.37  CO  +  2.29  H20  (volatile  partial  combustion) 

2.119  x  1011 

0 

2.027  x  108 

C1.37H4.5gO0.44  (volatile)  +  2.295  02  1.37  C02  +  2.29  H20  (volatile  complete  combustion) 

2.119  x  1011 

0 

2.199  x  1011 
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Fig.  2.  Meshed  geometry  (left),  closer  view  of  bottom  nozzles  (right). 


Fig.  3.  Grid  sensitivity  study:  temperature  (K)  distribution  along  height  (m)  for  five  grid  sizes. 


Table  4 

Operated  parameters  for  simulated  cases. 


Case  1st  Level  (upper  level,  AAMevel)  2nd  Level  (down  level,  BBMevel) 


name 

Fraction  of  total 
coal 

(%) 

Mass  flow  of 

coal 

(kg/s) 

Fraction  of  total 

oxygen 

(%) 

Mass  flow  of 

oxygen 

(kg/s) 

Fraction  of  total 
coal 

(%) 

Mass  flow  of 

coal 

(kg/s) 

%  Of  total 

oxygen 

(%) 

Fraction  of  total 

oxygen 

(kg/s) 

C3CL040 

30 

0.0006 

40 

0.0008 

70 

0.0014 

60 

0.0012 

C4CL040 

40 

0.0008 

40 

0.0008 

60 

0.0012 

60 

0.0012 

C5CL040 

50 

0.001 

40 

0.0008 

50 

0.001 

60 

0.0012 

C6CL040 

60 

0.0012 

40 

0.0008 

40 

0.0008 

60 

0.0012 

C7CL040 

70 

0.0014 

40 

0.0008 

30 

0.0006 

60 

0.0012 

C3CL050 

30 

0.0006 

50 

0.001 

70 

0.0014 

50 

0.001 

C4CL050 

40 

0.0008 

50 

0.001 

60 

0.0012 

50 

0.001 

C5CL050 

50 

0.001 

50 

0.001 

50 

0.001 

50 

0.001 

C60_050 

60 

0.0012 

50 

0.001 

40 

0.0008 

50 

0.001 

C70_050 

70 

0.0014 

50 

0.001 

30 

0.0006 

50 

0.001 

C30_060 

30 

0.0006 

60 

0.0012 

70 

0.0014 

40 

0.0008 

C40_060 

40 

0.0008 

60 

0.0012 

60 

0.0012 

40 

0.0008 

C50_060 

50 

0.001 

60 

0.0012 

50 

0.001 

40 

0.0008 

C60_060 

60 

0.0012 

60 

0.0012 

40 

0.0008 

40 

0.0008 

C70_060 

70 

0.0014 

60 

0.0012 

30 

0.0006 

40 

0.0008 

Note:  the  case  (C5CL060)  is  the  case  at  original  experimental  conditions  [37]  and  the  results  of  that  case  will  be  used  to  validate  the  model. 
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Fig.  4.  The  comparison  of  various  preliminary  simulation  results  with  experimental 
values  [37]  for  CO,  H2  and  C02  mol  fraction  in  syngas  (a)  at  O/C  ratio  =  0.9,  (b)  at  O/C 
ratio  =  1.0  and  (c)  at  O/C  ratio  =  1.1. 


2.4.  Boundary  conditions  and  calculation  methods 
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A  3D  computational  domain  has  developed  with  five  different 
grid  densities  to  analyze  the  grid  sensitivity.  Fig.  2  shows  meshed 
domain  along  with  the  closer  view  of  nozzles.  Temperature  is  the 
most  critical  parameter  in  gasification  so  it  was  selected  in  grid 
sensitivity  analysis.  Cold  flow  simulations  (no  reaction)  were  car¬ 
ried  out  with  all  five  different  grids  and  temperature  profile  along 
the  axis  for  all  the  grids  are  shown  in  Fig.  3.  The  grids  with  135,934 
and  256,223  tetrahedron  cells  shows  almost  same  temperature 
distribution  at  the  center  line  along  height  of  gasifier.  This  shows 
that  the  solution  is  independent  to  the  grid  beyond  135,934  grid 
size  hence  this  grid  has  selected  for  further  calculations. 


Fig.  5.  The  temperature  profile  for  various  preliminary  simulated  cases  and 
experimental  work  along  central  axis  of  gasifier. 


The  mass-flow  inlet  and  pressure-outlet  boundary  conditions 
were  used  for  all  input/output  streams.  Buoyancy  force  has  consid¬ 
ered  in  the  present  model.  Water  cooled  walls  were  assumed  to  be 
at  constant  temperature  at  800  K.  No-slip  state  (zero  velocity)  is 
applied  on  the  surfaces  of  walls.  Steady-state  simulations  were 
carried  out  with  an  implicit  pressure-correction  scheme 
(pressure-based  solver)  by  decoupling  energy  and  momentum 
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Fig.  6.  The  velocity  vectors  at  the  sectional  planes  of  the  gasifier  (Case  E)  with  a  close  view  of  AA'  and  BET  planes. 
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Fig.  7.  The  particles  residence  time  at  the  sectional  planes  of  the  gasifier  (Case  E) 
(for  clarity  only  200  particle  tracks  are  shown  with  particle  end  time  limited  at 
0.8  s). 

equations.  For  coupling  the  velocity  and  pressure,  SIMPLE  algo¬ 
rithm  was  followed.  Convective  terms  were  spatially  discredited 
by  second-order-upwind  scheme.  The  values  for  temperature 
dependent  properties  were  calculated  using  piecewise-polynomial 
equations  for  all  gas  and  solid  species.  Convergence  of  the  solution 
was  achieved  when  the  mass,  turbulent  kinetic  energy  and 
momentum  residuals  satisfied  at  1CT3  and  residuals  for  energy 
and  radiation  at  10~6.  Parallel  processing  was  used  for  computa¬ 
tion.  The  cold  flow  simulations  were  converged  first  in  all  the  cases 
then  the  reacting  flows  were  solved  by  activating  all  reactions 
along  with  the  injection  of  coal  with  ignition  temperature  of 
2000  K. 


3.  Results  and  discussion 

The  effects  of  coal  and  oxidant  distributions  for  two  stages  were 
investigated  through  various  simulations  by  varying  the  mass  frac¬ 
tion  (%)  of  total  coal  or  total  oxidant  at  both  injection  levels.  The 
overall  oxygen/coal  (O/C)  ratio  was  maintained  at  0.9,  1.0  and  1.1 
for  preliminary  validation  study  as  per  available  experimental  con¬ 
dition  [37].  The  rest  of  cases  were  simulated  through  validated 
model  at  constant  O/C  ratio  of  1.0.  The  names  of  simulation  cases 
referred  to  the  percentages  of  coal  and  oxygen  at  top  level  (AA'- 
level).  fable  4  describes  parameters  of  simulated  conditions  of  all 
cases. 

3.1.  Selection  of  best  reaction  mechanism  and  validation  of  the  model 

Fig.  4  shows  the  comparison  between  the  experimental  results 
[37]  and  various  preliminary  simulation  cases  (Case  A  to  F  from 
Table  2)  on  the  basis  of  mol%  of  CO,  H2  and  C02.  In  all  those  cases 
the  coal  or  oxygen  distribution  between  the  two  stages  is  kept  con¬ 
stant  as  per  experimental  conditions  (Coal  50%  (of  total)  at  AA' 
level  and  02  60%  (of  total)  at  AA'  level).  Three  O/C  ratios  0.9,  1.0 
and  1.1  were  taken  to  validate  the  model.  The  temperature  profiles 
along  the  axis  of  gasifier  for  all  simulated  cases  and  experimental 
observations  are  shown  in  Fig.  5  at  different  O/C  ratios. 

It  is  clear  from  Fig.  4  for  all  O/C  ratio  cases,  that  there  is  great 
impact  for  reaction-sets  on  the  overall  composition  of  syngas  com¬ 
position.  Case  A  follows  the  most  common  reaction-sets  i.e.  the 
partial  combustion  of  char  and  volatiles  and  the  involvement  of 
CO  combustion  in  the  bulk  gas  phase.  But  from  Fig.  4(a-c),  it  shows 
that  it  predicts  quite  less  CO  and  high  C02  as  compared  to  experi¬ 
mental  values  for  this  particular  type  of  gasifier  with  different 
arrangements  of  injecting  nozzles.  The  possible  cause  is  that  the 
well-mixed  pattern  of  particles  with  gas  due  to  impinging  and  tan¬ 
gential  injection  patters  of  fuel  and  oxygen.  The  velocity  vectors 
and  particles  are  shown  in  Figs.  6  and  7  that  gives  a  clear  idea  of 
flow  behavior  of  gas  inside  the  gasifier  along  with  the  250  colli¬ 
sions  of  particles.  Fig.  7  shows  the  particles’  residence  time  for 
200  tracks  only  with  limited  end  time  i.e.  0.8  s.  Fig.  7  confirms 
the  idea  that  there  is  a  greater  chance  to  have  the  particle  combus¬ 
tion  as  compared  to  CO  combustion  due  to  homogeneous  temper¬ 
ature  mixing  but  in  case  A  the  false  unrealistic  prediction  is  visible 
due  to  rapid  CO  combustion  and  increasing  the  temperature  of 
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Fig.  8.  The  mol%  for  CO,  H2  and  C02  with  variation  in  total  coal/oxygen  %  at  AA'  level. 
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Fig.  9.  The  char  conversion  with  variation  in  coal/oxygen  %  at  AA'  level. 


—■—40%  02 
—a—  50%  02 


Fig.  10.  Syngas  average  exit  temperature  (K). 


overall  gas  mix.  The  excess  temperature  rise  can  also  be  verified  in 
Fig.  5  for  same  case.  Rest  of  the  cases  (from  B  to  F)  shows  a  good 
agreement  of  H2  mol%  but  there  are  significant  differences  for  CO 


and  C02  mol%.  In  this  regard,  Case  E  shows  the  best  results  com¬ 
pared  with  experimental  values  for  prediction  of  syngas  composi¬ 
tion  and  temperature  at  different  O/C  conditions.  The  reaction 
mechanism  of  case  E  contains  the  complete  combustion  of  volatiles 
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and  char  and  ignoring  the  CO  combustion  reaction.  The  values  of 
CO  and  H2  mol%  for  this  case  are  in  good  agreement  and  show  less 
than  1%  error  while  the  prediction  of  C02  is  moderate  but  good  as 
compared  to  other  cases.  The  one  reason  for  not  predicting  the 
good  C02  is  the  assuming  no  formation  of  other  species  like  CH4, 
etc.  Temperature  prediction  through  case  E  is  also  in  satisfactory 
limits  (less  than  1%  error)  so  Case  E  has  been  validated  and  reaction 
mechanism  of  Case  E  has  selected  to  proceed  for  rest  of  studies. 

3.2.  Effects  of  coal/oxygen  distribution  on  syngas  composition 

15  Different  cases  were  simulated  at  O/C  ratio  =  1.0  by  varying 
coal  and  oxygen  distribution  at  both  the  injection  levels  as  per 
Table  4.  The  effects  on  the  mol%  of  CO,  C02  and  H2  by  varying 
the  total  coal  %  at  AA'  level  are  shown  in  Fig.  8.  As  per  Fig.  8(a), 
the  mol%  of  CO  shows  first  increasing  and  then  decreasing  trend 
by  increasing  coal%  at  AA'  level  with  40%  and  50%  oxygen  at  same 
level.  This  behavior  is  slightly  changed  with  60%  oxygen  case 
where  it  is  almost  of  increasing  order.  The  optimized  conditions 
for  maximum  CO  production  (52.59%)  was  found  at  50%  coal  and 
oxygen  at  AA'  level  (C50_050  case).  The  minimum  CO  mol%  was 
observed  50.123%  in  C30_060  case.  Exactly  inverse  trends  can  be 
seen  for  C02  mol%  with  respective  cases  (Fig.  8(b)).  A  plausible 
cause  for  these  trends  is  the  variation  of  oxygen/coal  ratio  at  local 
levels  and  the  rate  of  mixing  due  to  vortex  formation  with  the  tan¬ 
gential  nozzles  (either  at  AA'  level  or  BB'  level)  as  explained  by  Seo 
et  al.  [28].  The  maximum  and  minimum  for  C02  mol%  were 
observed  19.98%  and  18.39%  in  C30_060  and  C50_040  cases 
respectively. 


The  production  for  H2  was  found  in  the  range  of  27-28%  in  all 
the  cases  as  per  Fig.  8(c).  In  this  narrow  range  (27-28%)  it  shows 
a  bit  increasing  and  then  decreasing  trend  with  increasing  coal  at 
AA'  level  for  50%  oxygen;  a  decreasing  trend  with  40%  oxygen 
and  almost  increasing  behavior  with  60%  at  AA'  level.  The  reason 
for  this  behavior  of  H2  production  in  these  various  cases  is  the  non¬ 
availability  of  abundant  water  because  all  the  gasification  simula¬ 
tions  were  carried  out  without  steam.  Further  there  is  also  less 
moisture  (1.16%)  present  in  the  coal  as  per  its  proximate  analysis 
(refer  Table  1).  The  less  presence  of  water  has  a  great  influence 
on  the  water-shift  reaction  as  per  earlier  studies  [47,48]. 

3.3.  Effects  on  char  conversion 

The  coal/oxygen  variation  at  both  distributions  has  impact  on 
its  local  ratio  distribution  and  that  has  also  impact  on  the  overall 
char  conversion.  This  impact  is  shown  in  Fig.  9  for  various  simu¬ 
lated  cases.  In  general  an  increase  can  be  seen  in  conversion  by 
increasing  the  coal  at  top  level  (from  30%  to  50%)  with  all  oxygen 
cases.  But  conversion  decreases  by  increasing  coal  from  50%  to 
70%  at  top  level  (AA')  with  40%  oxygen  at  that  level.  The  other 
two  conditions  (50%  and  60%  oxygen  at  top)  show  a  continual 
increase  in  conversion  by  increasing  coal  from  50%  to  70%.  This 
behavior  is  due  to  greater  residence  time  of  the  coal  injected  at 
top-level  (AA'  level)  as  compared  to  the  coal  injected  at  bottom 
level  (BB'  level)  according  to  the  conclusions  of  earlier  studies 
[18,49].  So  the  higher  amount  of  coal  injected  at  AA'  level  with  suf¬ 
ficient  oxygen  has  enough  residence  time  to  reach  at  maximum 
conversion,  usually  more  than  99%.  The  minimum  char  conversion 
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Fig.  11.  Radial  temperature  profiles  for  various  simulation  cases  at  different  heights  of  gasifier. 
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i.e.,  95.45%  was  found  with  30%  coal  and  40%  oxygen  at  top  level 
case  (C30_040)  whereas  the  maximum  char  conversion  was  found 
99.7%  in  C60_050  and  C70_060  both  the  cases. 

3.4.  Effects  on  syngas  exit  temperature  and  maximum  inside 
temperature 

The  temperature  plays  an  important  role  in  the  gasification 
system.  The  variations  in  the  coal  and  oxygen  at  both  the  injection 
levels  actually  affect  local  oxygen/coal  ratio  as  discussed  in 


previous  sections  and  hence  the  defined  combustion/gasification 
reactions  occurred  with  different  rates.  Because  of  the  difference 
of  those  exothermic  and  endothermic  reaction  rates  there  is  an 
increase  or  decrease  in  the  overall  temperatures.  Fig.  10  shows 
an  overall  picture  for  the  syngas  exit  temperatures  for  simulated 
cases.  It  is  clear  that  in  all  the  cases  the  exit  temperature  for  syngas 
is  in  the  range  from  1250  I<  to  1450  K.  Li  et  al.  [34]  observed  that  in 
the  center  of  impinging  zone  there  is  an  increase  of  particle  cohe¬ 
sion  and  agglomeration  by  resulting  a  high  particle  concentration 
region.  The  rapid  deceleration  of  particles  near  the  center  of 
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Fig.  12.  Temperature  contours  for  sectional  planes  at  AA'  level  and  BB'  level. 
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Fig.  13.  The  particles  residence  time  at  the  sectional  planes  of  the  gasifier;  Cases:  C30_050,  C50_050  and  C70_050  (for  clarity  only  20  particle  tracks  are  shown  with  particle 
end  time  limited  at  0.2  s). 
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Fig.  14.  Turbulent  intensity  (%)  contours  for  selected  cases. 
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Fig.  15.  Net  heat  generated  due  to  reactions  (W)  for  selected  cases. 


impinging  zone  improves  the  performance  of  gasifier  in  terms  of 
heat  and  mass  transfer. 

The  radial  temperature  profiles  from  center  can  be  seen  in 
Fig.  1 1  for  three  coal  distribution  scenarios.  It  is  evident  that  there 
is  great  impact  for  the  local  temperature  distribution  due  to  varia¬ 
tion  in  coal  or  oxidant  between  two  stages.  Further  the  results  can 
be  verified  by  examining  the  temperature  profile  contours  of  top 
and  bottom  injection  sectional  planes  as  given  in  Fig.  12. 

It  is  apparent  from  those  temperature  contours  that  the  inside 
temperature  increases  by  increasing  the  amount  of  coal  at  any 
level  with  appropriate  amount  of  oxygen.  The  maximum  tempera¬ 
ture  was  observed  2027  K  at  AA'  level  with  70%  of  total  coal  and 
60%  of  total  oxygen  (Case  C70_060,  Fig.  12(c)). 


3.5.  Effects  of  coal  distribution  on  particle  trajectories 


The  particle  residence  time  is  key  parameters  for  the  conversion 
of  char  and  occurrence  of  different  reactions  with  their  specific 
rates,  as  explained  in  earlier  section.  The  effects  on  the  particle  res¬ 
idence  time  can  be  depicted  from  Fig.  13,  where  particles  trajecto¬ 
ries  are  shown  for  variation  of  coal  distribution  with  fixed  50% 
oxygen  distribution  at  both  the  levels.  For  clarity  the  trajectories 
are  limited  to  20  particle  streams  with  maximum  0.2  s.  The  varia¬ 
tion  in  coal  at  any  injection  level  with  fixed  oxygen  actually 
impacted  on  the  local  particle/gas  volume  ratio  and  hence  ulti¬ 
mately  impacted  on  the  particle  trajectories.  It  is  observed  that 
the  less  particle/gas  ratio  increases  the  particles  residence  time 
and  confirms  the  previous  studies  [44  . 
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Fig.  16.  Heat  consumption  in  different  ways. 


3.6.  Effects  on  turbulent  intensity 

Turbulent  intensity  has  predicted  by  solving  the  equations  of 
reliable  k-e  model  from  Eqs.  (5)-(9).  Due  to  impinging  and  tangen¬ 
tial  nozzles  at  AA'  and  BB'  levels,  the  major  turbulence  has 
observed  on  the  sectional  planes  on  these  levels.  The  effects  of  coal 
and  oxygen  distribution  on  Turbulent  Intensity  can  be  visualized 
from  Fig.  14.  The  turbulent  intensity  of  108-300%  has  observed 
at  the  main  reaction  zone  on  both  the  planes.  The  obvious  reason 
is  the  high  mixing  rate  at  these  sectional  planes  for  injectors. 

3.7.  Heat  generation  and  consumption  analysis 

The  overall  energy  balances  accounted  through  Eq.  (3).  The  pri¬ 
mary  source  for  generation  of  heat  energy  is  from  exothermic  reac¬ 
tions.  Then  the  heat  is  consumed  by  endothermic  reactions  and 
remaining  heat  is  either  absorbed  by  walls  which  are  assumed  to 
be  at  constant  temperature  of  800  I<  (water  cooled  system)  or 
taken  out  with  syngas.  A  fraction  of  heat  is  also  consumed  by  the 
solid  fuel  particles  to  raise  their  temperature  for  drying,  devolatil¬ 
ization  and  then  reactions.  The  net  heat  generated  through  reac¬ 
tions  has  observed  in  the  range  of  3000  W  and  shown  in  Fig.  15 
for  few  selected  cases.  The  maximum  heat  consumed  by  the  solid 
fuel  particles  and  observed  in  the  range  of  18,700  W  (Fig.  15).  Fleat 
lost  through  walls  is  found  in  the  range  of  9000  W  whereas  the 
heat  carried  away  by  syngas  is  found  in  2000-26,000  W  (Fig.  16). 

4.  Conclusion 

Double-stage  entrained  flow  gasifier  with  multi  opposite  burn¬ 
ers  for  lignite  coal  has  numerically  simulated  through  commercial 
CFD  software  in  present  study.  The  fuel  feeding  injectors  were 
impinging  and  tangential  in  nature  to  produce  slightly  vortex  flow 
inside  the  gasifier.  In  the  first  part  different  reaction  mechanisms 
were  simulated  as  per  experimental  conditions.  After  the  valida¬ 
tion  of  model  with  one  best  reaction  mechanism,  further  investiga¬ 
tions  were  made  to  study  the  effects  of  coal  and  oxygen 
distributions  between  the  two  feeding  injection  levels.  As  per  pre¬ 
liminary  simulations,  it  was  concluded  that: 

•  The  complete  combustion  of  char  and  volatiles  reaction  mecha¬ 
nism  can  predict  the  mol%  of  CO,  H2  and  C02  with  less  than  1% 
error. 

•  The  variation  of  oxygen  or  coal  at  any  injection  level  basically 
impacted  on  the  local  oxygen/coal  ratio  and  the  gas-solid  mix¬ 
ing  pattern.  This  impact  is  transferred  in  the  variation  of  char 
consumption  rates  with  combustion  and  gasification  reactions 
and  that  plays  a  key  role  in  the  variation  of  syngas  components 
and  temperature. 

•  A  little  variation  was  seen  in  H2  mol%  by  changing  coal  or  oxy¬ 
gen  at  both  the  levels  due  to  very  less  moisture  in  the  coal  along 
with  the  consideration  of  non-steam  gasification.  The  produc¬ 
tion  for  H2  was  found  in  the  range  of  27-28  mol%  in  all  the 
cases.  The  maximum  CO  mol%  was  found  52.59%  for  50%  coal 
and  50%  oxygen  at  upper  injection  level  (AA')  whereas  the  min¬ 
imum  CO  mol%  was  observed  22.37%  for  30%  coal  and  70%  oxy¬ 
gen  at  upper  level  (AA').  The  maximum  and  minimum  for  C02 
mol%  were  observed  47.86%  and  18.39%  in  30%  coal  with  70% 
oxygen  and  50%  coal  with  40%  oxygen  at  AA'  level  respectively. 

•  The  exit  temperature  for  syngas  was  found  in  the  range  from 
1250  K  to  1450  K.  The  maximum  temperature  was  observed 
2027  K  at  AA'  level  with  30%  of  total  coal  and  70%  of  total  oxy¬ 
gen  at  upper  injection  level. 

•  The  maximum  char  conversion  was  found  99.79%  with  coal  60% 
and  oxygen  50%  of  total  at  AA'  level.  The  minimum  char  conver¬ 
sion  was  observed  95.45%  at  30%  coal  with  40%  oxygen  at  AA'  level. 
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•  The  maximum  turbulent  intensity  has  observed  300%  at  the 
main  reaction  zone  on  both  injecting  planes  due  to  high  mixing 
rate  at  these  sectional  planes  for  injectors. 

•  In  general  with  oxygen/coal  above  or  equal  to  50%  of  total  at 
upper  injection  level  has  shown  an  optimized  performance. 

Overall  it  was  concluded  that  the  coal  and  oxygen  distribution 
has  great  effect  on  the  syngas  composition,  char  conversion  and 
exit  syngas  temperature.  The  temperature  and  kinetics  of  gasifica¬ 
tion  reactions  (reactions  of  char  to  C02  and  H20)  can  be  controlled 
with  the  optimized  coal  and  oxidant  distribution  between  the  two 
stages.  These  parameters  are  critical  to  the  overall  performance  of 
the  gasifier,  so  coal  and  oxygen  feedings  must  be  optimized 
between  the  two  stages  of  the  gasifier  to  get  the  optimized 
performance. 
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